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IV is  -shown- that  the  magnetic  moment  of  a  steel  ship,  whether 

permanent  or  induced,  can  conveniently  be  represented  by  an  array, 

*.  ,  ,  . 

inx three  dimensions,  of  magnetic  dipoles  at  the  centers  of  rectangu- 

•  #  ,  • 
lar  cells  which  the  ship  wholly  or  partly  occupies,  £he  magnetic  moment 

of -each  dipole  being  equal* to  that  of  the  portion  of  the  ship* s  volume 

which  -it  represents.  By-thiyartlflce  §  single  computation  of  magnetic 

field  components  of  a  magnetic  dipole  at  points  of  a  grid  centered  thereof, 

at*  can  be  used  to  find  all  three  magnetic  components  of  the  field  at 

any  . grid  point. dee -to  the  whole  ship  on  any  magnetic  heading  anywhere. 

The' agreement  with  experimental  data  should  be  good  at  points  much 

.  '  .f.«n  n.  .<■ 

closer  to  the  ship ,than  when  its  magnetic  moment  is  represented  by  a 

linear  distribution  of  magnetic  poles.  An  example  is' computed  in 
■  * 

some  detail. 
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'COMPUTATION  OF  SHIPS?  MAGNETIC  FIELDS 

•■•**•*,  •  >  •  • 

It  has  been  known  since  the' development  of  ;degatissingj  *durJ  ng 
World  War  II,  that  the  magnetic  field  of  a  ship,  at  points  where  the 

effects  of  its  local  inhomogeneities,  are  not  great;  is  nearly  .that  of  ■ 

,  '  *  •  ^  *  «• 

a  homogeneously  njagnetized  ellipsoid,  the  principal.  dimensions  of  which 

*  •  4  ,  ■  '  ’  ,  ‘ 

•  ,  ■  *  •  ♦  1 »’ 

are  nearly  the  principal  dimensions  of  the  ship  (length,  beam,  depth.),' 

and  which  is  centered  at  the  shipfs  centroid.  The  most  complete  : 

* 

statement  of  this  generalization  is  found  in  a  publication  (1)  of  the  • 
Bureau  of  Ordnance  hereinafter  referred  to  as  OD  £050.  ‘ 


.  In  OD  $050,  (pp..  432*463,  in  Part  IV)  a  methyl  u-lr.g  spheroidal 

coordinates,  applicable  where  the  more  general  ellipsoid  degenerates" 

*  •  •  •  * 

into. a  spheroid,  is  explained,  and  is  applied  to  the  USS  OVERTON  DD  239, 

to  show  how  magnetic  measurements  made  on  this  ship  in  1949  can  be  repro#* 

« 

.duced  by  computation.  It  is  correctly  stated  that  such  computations  are 
inconvenient,  and  are  too  time-consuming  for  gneral  use. 


In  another  place  in  OD  $050  (p.  69,  in  Part  I)  an  empirical  formula 

is  given  for  finding  a  vector  proportional  to  the  magnetic  moment  of 

0 

a  ship  heading  in  any  direction  in  the  local  magnetic  field  of  the  earth. 
This  magnetic  moment, or  M,  can  safely  be  used  to  compute  the  magnetic 
field  of  the  ship  at  distances  great  with  respect  to  its  dimensions. 

It  is  the  purpose  of  this  note  to  point  out  that , it  is  possible, 

•  • 

and  not  too  time-consuming  if  high  speed  digital  computers  are  available, 
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to  combine  the  two  idsas  just  presented,  and  to  compute. relatively 
near  fields  of  a  ship,  considered  as  a  uniformly  magnetized  ellipsoid,' 
by  adding  up  the  separately  computed  field  components  of  suitable  . 
transverse  slices  of  the  ellipsoid.  A  trial  computation  for  the  OVERTON, 
done  by  "hand",  using  a  desk  type  ■  electrically  driven  digital  machine.,. 
ha3  beejfi  carried  through  and  will  be  presented  below.  The  use  of  the 
new  method  for  predicting  the  expected  interaction  'of  steel  ships  of 
known  dimensions  with  magnetic  mines  of  known,  or  assumed,  character¬ 
istics,'  may  turn  out  to  be  justified  in  solving  complicated  mine  war¬ 
fare  problems. 

.  .  *  ¥  * 

The  formula  for  M,  given  in  OD  £050,  is  , 

.  V  M  s  0.013X(1*- p)i  ♦  0.003Tj  *0.005Zk  (1) 

* 

•  Here’ X,Y,Z  are  components  of  the  earth’s  local  magnetic  field 
along  the  unit  vectors  fixed  in  the  ship,  pointing  forward,  £  .. 

pointing  to  starboard,  k,  pointing  downward.  Thus  X  and  Y  depend  upon 
the  ship’s  heading,  but  Z  does  not.  The  number  p  is  stated  to  be  "the 
ratio  of  permanent  to  induced  longitudinal  magnetization",  but  study 
of  the  formula  shows  that  this  is  not  a  possible  definition  for  p 
since  the  formula  would  then  say  that  when  X  becomes  zero  th*e  "permanent" 
magnetization  has  no  effect  on  JJ,  which  is  absurd.  The  mistake  must 
have  been  committed  in  simplifying  a  general  formula  for  use  with  con¬ 
ventional  north  and  south  headings  in  degaussing  work.  The  general 
formula,,  which  should  have  been  given,  is  * 
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*  ,■  * 

'  M--  0,013  (X+ P)i  0.003  Y£+0.005Zk,  (2) 

*  t  «  • 

* 

in  which  P  stands  for  a  magnetic  field,  constant  in  magnitude  and  in  ■ 

direction  along  which  would  account  for  the  permanent  longitudinal 
.  •  • 

magnetization  ad  X  accounts  for  the  induced  longitudinal  magnetization. 

‘A  more  general  expression  of  Mg,  the  magnetic 'moment  of  the  ship, m 

*  m  *» 

i3  *  •  . 


Mg  =  Kx  (x.+  P)  ltKy  Yj4-Kz  Zk  (3) 

where  Kx,  Ky,  Kz  are  factors,  all  of  the  same  physical  dimensions 
(magnetic  moment  divided  by  magnetic  field  intensity).  Since  remote 
field  changes,  radial  and  circumferential,  tan  be  computed  from 

A  H*  s  2Mg  cos  e/,,3  (4)  ..  . 

4  H«  .  Mg  sin  ©4 3  (5) 

it  is  clear  that,  in  Gaussian  units,  here  used,  the  physical  dimensions 
of  magnetic  moment  divided  by  magnetic  field  intensity  are  those  of  ♦ 
volume,  or  This  weans  that  here,  as  in  OD  3050,  lengths  can  be 

expressed  in  reduced  form,  that  is  as  numerical  multiples  of  any  con¬ 
venient  length,  in  both  parts  of  the  computation  leading  to  magnetic 
field  values.  It  also  states  the  well-known  fact  that  scale  models  can 
be  used  in  planning  the  degaussing  of  ships.  It  is  also  clea.  that* all 
magnetic  fields  can  be  expressed  in  oersteds  or  in  mill'i-oersteds,  or, 
as  is  customary  in  degaussing  circles,  in  gausses  or  milli-gausses. 

(The  use. of  units  of  magnetic  induction  is  tolerable  because  ships  and 


CONFIDENTIAL 


-k- 


CONFIDENTIAL* 

r  .  i 

% 

mines  interact  in  njedia  which  cannot  be  appreciably  magnetized.) 

Compuatation  by  ship  slices  can  be  much  shortened  by  choosing  for 
computation  points  in  the  field  which  are  separated,  along  the  direction 
of  the  keel,  by  intervals  equal  to  the  thickness  of  each  slice.  The 
numerical  factors  which  appear  in  computing  the  field  compononts  due  to 
one  slice  at  a  longitudinal  row  of  points  will  reappear  for  every  other 

slice  at  other  points  of  the  row  and, .so  need,  be  computed  only  once. 

•  , 

The  expansion  of  the  vectorial  formulae  for  field  component  of  the  n^ 
slice,  with  its  resultant  magnetic  moment  FnMn  at. its  centroid,  yields 
the?' following  formulae  for  a  well-depermed  ship,  for  which,  that  is, 

-  4  '  . 

P-0*. 

v 

...  * 

i  *, 

AHxn  r  Fn  (Kx^X+Ky*'Y*K2/rf'2i)  (6) 

4Hyn  =  Fn  (K^'X+Ky^Y+it^'z)  (7) 

^Hzn  :  Fn  (Kx^^Ky^'Z^Kz^Z)  (3) 

•  • 

Here  Fn  is  the  fraction  of  the  volume  of  the  ellipsoid  in  the  slice  and 
the  factors  of,  /$ ,  ){,  fat*  are  numerical  factors  expressed  in  the 

reduced  coordinates  V>,  ^  of  the  point  chosen  with  respect  to  the 
centroid  of  the  slice,  and  of  the  reduced  radius  from  that  centroid. 
Explicitly 

f  *  (9) 

-jn  (io (ii)<*-~t-~t  (12) 
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-  3  \  ^  /  3  r  / ,  ,  *  v* •  5  %  *-> 

*  '  ~ZZrr  {1*‘/S-m  ^  .(u)  <*>  ~^T~  ■■ 


(15) 


The  relative  magnitudes  of  Kx,  Ky,  Kz  'here  assumed,  were  consistent 
with  CD  3050,  that  is 


.  Kx  *  Ky  :  Kz  r’O. 013:0. 003:0.005  (16)  - 

The-  common  multiplier  of  the  numbers  on  the  'right  hand  .side  of  this  •  ’ 
proportion,  to  give  Kx,  Ky,  Kz  is  (  J? / 2)3  whehe.^'is  the  length ‘of  the 
ship  (not  the  length  of  the  substituted  ellipsoid)  For  the  OVERTON 
of  length  314  ft  an'd  for  the  chosen  unit  of  length,.  40'  ft;  it  follows 
that  the  reduced  length  is  A  2  7*35*  The • numerical  factor  used  in 
computations  is  then  (<A/2)^  r  60,467* 


The  unit  length,  40  ft,  wa3  chosen  so  that  computations  for  points 
on  a  plane  40  ft  below  the  ship’s  centroid  (as  in  OD  3050)'  would  be 
especially  easy.  It  was  thought  better,  however,  to  use  slices  less 
than  40  ft  thick,  and  compute  fields  for  points  le3s  than  40  ft  apart 

•1 

on  lines  parallel  to  the  ship’s  course.  It  can  now  be  concluded  that 
20  ft  would  have  been  about  right  as  the  slice  thickness  and  grid  inter¬ 
val,  but,  in  fact,  10  ft  was  chosen.  This  made  the  central  slices  of 
the  OVERTON  considerably  too  thin  to  be  well  represented  by  point  mag¬ 
netic  dipoles  at  their  centroids,  and  this,  in  turn,  makes  all  computed 
fields  at  f  :  1  too  samll  by  possibly  3 . per  cent  at  p&ak  values,  posi- 
tive  or  negative.  Less  than  peak  values  arc  proportionately,'  as  well  as 
absolutely,  closer  to  the  best  values  obtainable  by  the  slice  treatment. 

The  three  numbers  given  in  Eq.  (1)  and  copied  in  Eq.  (2),  are 
* 
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different  for  three  reasons.  Most  fundamentally,  the  3teel  in  a  ship 
is  not  isotropically  arranged,  there  being  different  average  areas  of 
steel  in  random, plane  sections  perpendicular  to  the  three  vectors  ^i,  ^ 

k.  Secondly,  the  resultant  magnetic  moments  induced  by  equal  axial 

*>  * 

fields  in  a  magnetically  isotropic  ellipsoid  would  still  diffor  along 

different  axes,  being  greater  in  the  directions  of  longer  principal 

axes  because  the  jelf-demagneti 2ing  fields  are  leas  for  then*  directions 

*  • 

Thirdly,  the  effect  of  permanent  magnetization  is  included  only  in  the 

third  term.  The  numbers  given  in  OD  3050  are  only  appropriate  for  ships 
of  conventional  shape  and  structure.  Any  radical  departure  from  the 
type  of  hull  on  which  these  numbers  were  based  would  change  their  rela¬ 
tive  magnitudes.  Ail  this  means  is  that  more  recent  experience  must  be 
utilized  before  any  extensive  computations  can  be  based  upon  Eq.  (2), 
and-ttiat  different  multipliers  Kx,  Ky.  Kz  can  be  expected  for  ships  of 
peculiar  characteristics. 

The  spheroid  found  in  CD  3050  to  be  best  suited  to  calculations  of 

s 

OVERTON  field  values  by  spheroidal  harmonic  analysis  was  also  assumed 
here.  It  has  foci  352  ft  apart  and  a  major  axis  353.2  ft  long.  It  was 
divided  into  35  slices,  the  central  slice  centered  at  x  *  0,  y  =  0,  z  * 
and  the  extreme  slices  centered  at  x  z  *  170  ft,  y  •  0,  z  -  0.  (Thi3 
ignores  a  very  small  volume  at  each  end.  )  Neglecting  the  curvature  of 
the  slice  edge3  the  relative  volumes  were  taken  as  proportional  to  the 
areas  of  their  transverse  mid-sretions.  The  resulting  values  of  **v  are 

JT 

given  in  Table  1,  and  have  been  adjusted  to  make  ^  yf  si* 
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Table  I 

Fractions,  ,  of  total  magnetic  moment  of  U3S  OVERTON  DD  239  ellipsoid 


model,  contained 

in  transverge  slice, w , 

10  ft  thick. 

Yl 

} 

*  ■  •  ■  a 

1 

35 

-4.25 

0.003114 

2 

34 

-4.00 

0.007607 

*  3 

33 

-3.75 

0.011327 

4 

32 

-3.50 

0%015775 

5 

31 

-3.25 

0.019450 

6 

30 

-3.00. 

0.022354 

7 

29 

-2.75 

0.025935 

3. 

23 

-2.50 

0.023344 

9 

27 

-2.25 

0.031430 

10 

26 

-2.00 

0.033745 

11 

25 

-1.75 

0.03573 7 

12 

24 

-1.50 

0.037556 

13 

23 

-1.25 

4 

0.039054 

14 

22 

-1.00 

0.040279 

15 

21 

^  .75 

0.041232 

16 

20 

-  .50 

0.041912 

17 

19 

2  .25 

0.042321 

13 

0 

0.042457 

# 
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Computation  is  straight-forward  for  any  point  (  $  ,  n  ,  ^  )  in  the 

field,  being  the  summation  for  each  of  the  nine  terms  in  Eqs.  (6),  (7), 

[&)  of  35  products  of  4  by  appropriate  *  , - ^'factors.  The 

values  of^  t/At  )f  v'  are  tabulated  for  a  few  fixed  values  of  ri 

(0,  -0.5,  -3,  -2,  -3)#  for  one  value  of  /(l),  and  for  a  range  in 

values  of  f (0,  *0.25,  -0.50, -  *9.75,  *10.00)  sufficiently  great  to 

include’ points  well  ahead  of  the  ship  su.d  well  behind  it.  Th«  chosen 

values  of./  andbpermitted  mapping  the  field  of  the  OVERTON  in  a  . 
t  * 

rectangle  460  ft  long  and  240  ft  wide,  or  to  points  73  ft.  ahead  and 
behind,  to  points  about  105  ft  off  the  beam.  The  assumed  values  of 
the.ear*£hVs  magnetic  field  were  H  (*(*a  *  /a)i)  0.200  gauss,  H(»  Z  ) 
0.49,5  cnt>3S,  as  at  the  point  where  the  OVERTON  was  studied..  All  H 
terms  directly  computed  (before  summing  the  35  items)  we-e  recorded  to 
10“^  gauss. 

0 

If  a  grid  of  points  at  the  corners  of  square  meshes  had  been  com¬ 
pletely'  filled  in,  by  selecting  the  range  and  intervals  in  Vi  to  match 
those  in  £ ,  the  computation  would  not  have  been  as  much  greater  as  a 
mere  count  of  added  grid  points  would  suggest,  since  symmetry  reduces 
the  task  considerably. 

Figs.  1-3,  4-6,  7-9  show  ,4  Hy  t  A  at  the  selected  grid 
points  for  the  OVERTON  headed  0°,  45°  and  90°  (magnetic)  where  the 
earth's  field  has  a  horizontal  component  of  0.209  gauss,  a  dcwn-point- 
ing  vertical  component  of  0.495  gauss. 

Since  there  are  mines  which  attend  only  to  changes  in  the  total 
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HT  r  (x2*y2*z2)*  «  (H2Vh2)*  (17) 

•  • 

.  • 

magnetic  field  it  seemed  desirable  to  compute  the  magnitude  of  the 

*  -  r<HT{!3) 

resultant  of  A  Hx,£Hy,.AH2  (at  each  point  for  which  these  components  • 
.had  been  computed).  This  was  done,  and  Figs.  10-12  present  the  results 
graphically  for  the  same  three  headings. 

e> . 

Results  for  other  cardinal  and  intercardinal  headings  can  easily 
be  derived  by  symmetry  from  those  given  in  the  Figures.  Interpolation 
for  other  headings  and  for  other  values  of  n ,  may  also  be  feasible. 

The  behavior  of  a  mine  attending  to  the  total  magnetic  field,  or 
to  changes  in  its  vertical  component,  when  a  ship  passes,  can  be  pre¬ 
dicted  from  curves  like  those  already  given.  Plotting  field  level 
‘lines  on  a  horizontal  projection  of  the  ship  and  its  neighborhood  to 
get  actuation  widths  for  mines  of  known  characteristics  is  a  well-known’ 
method  used  in  such  prediction.  Mines  attending  to  changes  in  the 
horizontal  component  need  a  further  step  in  computation,  explained 
below. 

A  horizontal-component  induction  type  mine  can  only  attend  to 
changes  in  Hc^the  horizontal  component  which  lies  along  the  axis  of  the 
core  of  its  induction  coil.  The  value  of  this  change  4H$,  at  any 
grid  point,  iadH^s  A  Hx  A  Hy  /a**,  $  when  ^  is  angle  between  ship 

heading  and  mine  axis.  It  can  be  assumed,  unless  the  contrary  can  be 
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proved,  that  all  orientations  of  this  axis  in  the  horizontal  plane  of  the 

sea  bottom  are  equally  probable.  If  this  be  admitted,  the  chance  that 

such  a  mine  will  be  fired  when  passed  by  a  ship  on  a  given  heading 

depends  upon  the  angle,  between  the  core  axis  and  the  magnetic 
* 

meridian,  as  well  as  u>on  the  ship»s  magnetic  heading,©. 

The  computation  of  the  changes  in  at.  grid  points*  for  which  A Hx 
andAHy  have  already  been  computed  is  again  straightforward.  It  has 
been  carried  out  for  the  came  •values  (0°,  45 90°)  for  a  series  of  # 
values  0°,  30°,  60°,  90°.  What  is  more  interesting  than  the  results 
for  any  specific  <)>,  for  a  given  ©  ,  is  the  mean  value  for  random 

orientations  01  the  mine  axes,  here  approximated  by 

=  (4  H0o ♦  2  A  H30o  * 2  A  H6qo  +A  H9Qo)/6 

« 

Values  ofAH^  for  9  a  0°,  450,  900  are  preoented  in  Figs.  13-1*5,  The 
contour  map  technique  may  again  be  used  to  get  most  probable  results  and 
to  estimate  their  variances  in  different  parts  of  real  channels. 

If  it  turns  out  that  single-component  mines  do  not  lie  close  to 
horizontal,  but  are  predominantly  at  an  appreciable  angle  to  the  hor¬ 
izontal  plane  when  waiting  for  ships,  the  analysis  here  used  will  have 
to  be  replaced  by  a  similar,  but  less  simple,  analysis. 

Another  method  of  handling  the  dependence  of  actuation  probability 
upon  the  orientation  of  horizontal  component  mines  was  developed  (l>) 

(3)  U)  in  connection  with  magnetic  mine  sweeping.  It  uses  vector 
diagrams  of  horizontal  field  increments  at  equal  time  intervals  at 
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fixed  points  on  the  botcora  as  the  field- producing  agent  passes,  and  a 

I 

transparent  rotatable  overlay  showing  the  assumed  orientation  of  the 
mine  axis  and  required  actuating  field  minima  along  that  aids,.  This 
procedure  is  time  consuming  as  first  proposed,  and  is  not  well  adapted 
to  machine  computing  even  in  its  latest  form,  Since  no  assumptions 
.were  made  in  the  present  paper  regarding  sensitivities  dr  other  operating 
characteristics  of  mines  this  method  was  here  inapplicable. 

Values  of  A H2,  for  0  s  0°,  for  various  values  of  x  and  for  y  3  0, 

y  z  20  ft,  are  replotted  in  Fig*  16  in  order  to  allow  comparison  with 

OVERTON  data  given  in  OD  8050  (on  pp .  453  and  454)  from  measurements  at 

the  Wolf  Trap  degaussing  station  near  Norfolk,  Va  These  data  being 

presented  in  the  form  of  contour  maps  for  0  ;  0°  and  0  *  180°,  it  was 

necessary  to  interpolate  between  contours  to  get  curves  of  AH  -r^-x 

z 

for  Fig.  16.  The  depth  to  magnetometers  was  not  exactly  40  ft,  measured 
from  the  ship  centroid  (assumed  to  be  2  ft  above  the  waterline)  but 
z  was  nearly  enough  40  ft  to  make  comparison  interesting. 

It  is  seen  that  the  course  of  values  is  very  similar  in  both, 
cases,  but  that  computed  values  are  smaller  at  the  maximum,  an*  show 
more  dissymmetry  about  the  midship  section  (x  -  0),  so  as  to  have  weak 
negative  values  near  the  stern.  These  differences  mean  that  in  the 
OVERTON  the  vertical  magnetization  is  greater,  and  the  longitudinal 
magnetization  is  less,  than  is  predicted  by  the  factors  given  for  ships 
in  general  early  in  OD  8050.  It  is  suggested  that  this  is  due  to  the 
fact  that  a  destroyer  has  more  vertical  steel  in  proportion  than  vesselo 
of  greater  beam  and  larger  compartments  in  which  decks,  bottom,  and  side 
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plying  are  relatively  heavier.  This  emphasizes  the  need  for  establish- 
in*  relative  K  -values  Eq.  (16)  for  ships  of  different  types  if*  cal¬ 
culated  magnetic  fields  of  ship3  are  to  be  relied  upon. 

Sirce  the  computations  here  reported  were  completed  I  have  found 

a  note  by  It.  K.  Reber  (5)  giving  a  detailed  analysis  of  the  magnetic 

# 

field  of  a  barge  110  feet  long,  30  feet  wide,  and  24  feet,  deep.  The 
width  and  depth  were  therefore  almost  the  same  as  for  the  OVERTON,  the 
length  only  one  third  as  great.  Reber  did  not  attempt  to  compute  the 
field  of  the  barge  at  nearby  points,  except  for  the  vertical  field  com- 
ponent  at.  a  point  on  the  bottom  under  the  center  of  the  vessel,.  He  says, 
inf:,. fact:  "Accurate  theoretical  calculations  of  the  horizontal  component 
fields  due  to  longitudinal  and  transverse  magnetization  are  hardly 

f' 

•  possible. "  He  does  estimate  the  total  longitudinal  magnetization  for  a 
no^th  or  3outh  heading,,  and  the  total  transverse  magi vtizat ion  for  an 
east  or  west  heading,  From  these  it  is  easy  to  find  relative  values  of 
Kx,  Ky,  Kz  in  the  notation  of  the  present  paper  In  order  to  make  com¬ 
parison  easier  I  have  arbitrarily  assumed  Kz  -  0.005  in  the  following 
equation  derived  from  his  findings. 

Kx  :  Ky  •  Kz  :  0.0054  :  0. 00.37  :  0.005  (20) 

It  is  seen  that  the  relative  values  of  Ky  and  K7  agree  excellently  with 
corresponding  numbers  in  Eq.  (16),  but  that  Kx  is  only  about  two  fifths 
of  the  earlier  given  value.  This  is  clearly  due  to  the  relative  short- 
. ness  of  the  barge  as  compared  with  typical  cargo  ships. 
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The  piaee  where  the  magnetic  measurements  uSed  *by  Reber  were  made 
was  very  near  where  fh-»  earlier  measurements  >n  the  OVERTON  were. made- 
Reber  assumed  fer  the  earth’s  magnetic  field  components  the  val'ue3 
Hv  ;  0.^  gauss.  0  1*  gantr .  whereas  the  values  used  in  tho  present 

paper  were  Hv  r  0-495  gauss,  i.  0,200  gauss.  If  the  present  cal¬ 
culations  were  *•:•  be  corrected  so  as  to  agree  with  Reber* s  assumptions 
M„  would  be  ftbmii  per  cent  greater,  irni  iA  f.  v  .  nor*  n  heading, 
would  be  10  per  cent,  less  Both  of  these  changes  would  improve  the  fit 
between  observed  and  calculated  curves  for  AH,  lr  Fig  16 

My  thanks  are  due  to  Mrs.  C  K  Bockelmnn  and  Mr3.  F.  H.  .Stillinger, 
Jr.  for  computation  f  all  field  component  values  used  in  the  Figures, 
except-  those  presented  as  "observed"  values  of  A  H?  in  Fig  16,  which 
I  derived  from  small  scale  vertical  field  contour  maps  In  Figures  15*4-5 
and  15.46  in  CD  3051),  and  which  are  much  less  precise 
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